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The effective-surface approximation is extended taking into account derivatives of the symmetry- 
energy density per particle with respect to the mean particle density. The isoscalar and isovector 
particle densities in this extended effective-surface approximation are derived. The improved expres¬ 
sions of the surface symmetry energy, in particular, its surface tension coefficients in the sharp-edged 
proton-neutron asymmetric nuclei take into account important gradient terms of the energy density 
functional. For most Skyrme forces the surface symmetry-energy constants and the corresponding 
neutron skins and isovector stiffnesses are calculated as functions of the Swiatecki derivative of the 
nongradient term of the symmetry-energy density per particle with respect to the isoscalar den¬ 
sity. Using the analytical isovector surface-energy constants in the framework of the Fermi-liquid 
droplet model we find energies and sum rules of the isovector giant-dipole resonance structure in 
a reasonable agreement with the experimental data, and they are compared with other theoretical 
approaches. 
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I. INTRODUCTION 

Explicit and accurate analytical expressions for the 
particle density distributions within the nuclear effective- 
surface (ES) approximation were obtained in Refs. [1-3]. 
They take advantage of the saturation properties of nu¬ 
clear matter in the narrow diffuse-edge region in finite 
heavy nuclei. The ES is defined as the location of points 
with a maximal density gradient. An orthogonal coor¬ 
dinate system related locally to the ES is specified by 
the distance ^ of a given point from the ES and the tan¬ 
gent coordinate 77 parallel to the ES. Using the nuclear 
energy-density functional theory, one can simplify the 
variational condition derived from minimization of the 
nuclear energy at some fixed integrals of motion in the 

r] coordinates within the leptodermous approximation. 
In particular, in the extended Thomas-Fermi (ETF) ap¬ 
proach [4], this can be done in sufficiently heavy nuclei 
for any fixed deformation using the expansion in a small 
parameter a/R ^ <C 1 where a is of the order of 

the diffuse edge thickness of the nucleus, R is the mean 
curvature radius of the ES, and A is the number of nu¬ 
cleons. The accuracy of the ES approximation in the 
ETF approach was checked [3] without spin-orbit (SO) 
and asymmetry terms by comparing results with those 
of Hartree-Fock (HE) and other ETF models for some 
Skyrme forces. The ES approach [3] was also extended 
by taking into account the SO and asymmetry effects 
[5-7]. 

Solutions for the isoscalar and isovector particle den¬ 


sities and energies in the ES (leptodermous) approxima¬ 
tion of the ETF approach were applied to analytical cal¬ 
culations of the surface symmetry energy, the neutron 
skin and isovector stiffness coefficient in the leading or¬ 
der of the parameter a/R [7]. Our results are compared 
with older investigations [8-11] within the liquid droplet 
model (LDM) and with more recent works [12-27]. 

The splitting of the isovector giant-dipole resonances 
into the main and satellite modes [28] was obtained as a 
function of the isovector surface-energy constant within 
the Fermi-liquid droplet (ELD) model [29, 30] in the 
ES approach. The analytical expressions for the sur¬ 
face symmetry-energy constants have been tested by en¬ 
ergies and sum rules of the isovector dipole resonances 
(IVDR) within the ELD model [31] for some Skyrme 
forces neglecting derivatives of the nongradient terms in 
the symmetry-energy density per particle with respect 
to the mean particle density. The so called pygmy dipole 
resonances (PDR) below the main IVDR peak as a differ¬ 
ent phenomenon which might not be the result of split¬ 
ting of the IVDR were intensively discussed in the liter¬ 
ature [12-14, 22-25, 32-35]. They might have a different 
nature and actually are not related to each other. In the 
present work, we shall extend the variational effective- 
surface method accounting for the derivatives introduced 
by Swiatecki and Myers within the LDM [8] and apply 
it to the IVDR splitting. Some preliminary results were 
reported in [36]. 

In Sec. II, we give an outlook of the basic points of the 
ES approximation within the density functional theory. 
The main results for the isoscalar and isovector particle 
densities are presented in Sec. Ill emphasizing derivatives 
of the symmetry energy density per particle. Section IV 
is devoted to analytical derivations of the symmetry en¬ 
ergy in terms of the surface energy coefficient, the neu- 
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tron skin thickness and the isovector stiffness including 
these derivatives. Sections V and VI are devoted to the 
collective dynamical description of the IVDR structure 
in terms of the response functions and transition densi¬ 
ties. Discussions of the results are given in Sec. VII and 
summarized in Sec. VIII. Some details of our calculations 
are presented in Appendixes A and B. 


II. SYMMETRY ENERGY AND PARTICLE 
DENSITIES 

We start with the nuclear energy E as a functional 
of the isoscalar (p_|_) and isovector (/?_) densities p± = 
Pn ± Pp in the local density approach [4, 37-44]: 

E = J dr p+S (p+, P-), (1) 

where £ {p+, p-) is the energy density per particle, 

^ (p+jP-) = —by + JI"^ + e+{p+) -I- S- {p+, P-) + 

+ {C+/P+ + V+) {Vp+f + (C_/p+ + D_) {Vp. f . (2) 

Here, by « 16 MeV is the separation energy of a particle, 
J « 30 MeV is the main volume symmetry-energy con¬ 
stant of infinite nuclear matter, and I = {N — Z)/A 
is the asymmetry parameter; N = J dr/9„(r) and 
Z = J drpp{r) are the neutron and proton numbers, and 
A = N + Z. These constants determine the first two 
terms of the volume energy. The last four terms are sur¬ 
face terms. The first two terms are independent of the 
gradients of the particle densities and the last two depend 
on these gradients. For the first surface term independent 
of the gradients, e+, one obtains 

e+ip+) = § e+ [e{w+)] , (3) 

where K+ « 215 — 245 MeV (see Table I) is the isoscalar 
in-compressibility modulus of symmetric nuclear matter, 
is the dimensionless isoscalar particle density, w+ = 
p+IV and 

e+ [e (w+)] = -f [5sym(e) - J] /K+ (4) 

with 

p-p+ 1-W+ 

e is the small parameter in the expansion, 

‘5sym(e) = J — Le -\—..., (6) 

around the particle density of infinite nuclear matter 
p = 3/47rrQ « 0.16 fm“^ and Tq is the commonly ac¬ 
cepted constant in the A^^^ dependence of a mean radius. 
Several other quantities, which were introduced by Myers 
and Swiatecki [8], will be explained below. The deriva¬ 
tive corrections of Ssymie) in Eq. (6) were neglected in 


our previous calculations [7]. The next isovector surface 
term e_(p+,p_) can be defined through the same func¬ 
tion iSsym(e) in Eq. (6): 

£_(p+,p_)=5.y„,(e) -J/2. (7) 


For the first and second derivatives of Ssym{f-) with re¬ 
spect to e one can take in Eq. (6) the values L « 20^ 120 
MeV and, even less known, K_ [19, 27, 45]. The con¬ 
stants C± and 'D± in Eq. (2) are defined by the parame¬ 
ters of the Skyrme forces [4, 37, 39, 41, 43], 


C+ = 
c_ = 


If 25 5 

12 r ■ ~ 3^"^" 


'48 


1 + - 


h 

36 


1 19 

1 + ^- 2^2 


( 8 ) 


The isoscalar SO gradient terms in (2) are defined with 
a constant: D+ = —9mlFQ/16?i^, where Wo « 100-130 
MeV-fm^ and m is the nucleon mass. The constant D_ 
is usually relatively small and will be neglected below for 
simplicity. Equation (2) can be applied in a semiclassical 
approximation for a realistic Skyrme force [37, 39-42], 
in particular by neglecting higher h corrections in the 
ETF kinetic energy [2-4] and also Coulomb terms. All 
of them easily were taken into account [1, 5] neglecting 
relatively small Coulomb exchange terms. Such exchange 
terms can be calculated numerically in extended Slater 
approximations [46]. 

The energy density per particle in Eq. (2) contains 
the first two volume terms and surface components in¬ 
cluding the new L and K- derivative corrections £_ (in 
contrast to Ref. [7]) and also the isoscalar and isovector 
density gradients. Both are important for finite nuclear 
systems. These gradient terms together with the other 
surface components in the energy density within the ES 
approximation are responsible for the surface tension in 
finite nuclei. 

As usual, we minimize the energy E under the con¬ 
straints of fixed particle number A = / dr (r) and neu¬ 
tron excess N — Z = J dr p_ (r) using the Lagrange mul¬ 
tipliers A+ and A_, the isoscalar and isovector chemical- 
potential surface corrections (see Appendix A). Taking 
also into account additional deformation constraints (like 
the quadrupole moment), our approach can be applied 
for any deformation parameter of the nuclear surface, if 
its diffuseness a is small with respect to the curvature ra¬ 
dius R. Approximate analytical expressions of the bind¬ 
ing energy will be obtained at least up to order 
To satisfy the condition of particle number conservation 
with the required accuracy we account for relatively small 
surface corrections (oc a/R ~ in first order) to the 

leading terms in the Lagrange multipliers [2, 3, 5, 7] (see 
Appendix B). We take into account explicitly the dif¬ 
fuseness of the particle density distributions. Solutions 
of the variational Lagrange equations can be derived ana¬ 
lytically for the isoscalar and isovector surface tension co¬ 
efficients (energy constants), instead of the phenomeno¬ 
logical constants of the standard LDM [8] (the neutron 
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and proton particle densities were considered earlier to 
be distributions with a strictly sharp edge). 


III. EXTENDED ISOSCALAR AND 
ISOVECTOR DENSITIES 


For the isoscalar particle density, w = p+ jji, one has up 
to the leading terms in the leptodermous parameter ajR 
the usual first-order differential Lagrange equation [3, 5, 
7]. Integrating this equation, one finds the solution: 


x = - dy 


I i + Py 

y e+[e(y)] ’ 


x = 

a 


(9) 


for X < x{w = 0 ) and ru = 0 for a; > x{w = 0 ), where 
x{w = 0 ) is the turning point. j3 = 'D+'p/C-[. is the di¬ 
mensionless SO parameter, see Eq. (4) for e+[e{y)] (for 
convenience we often omit the lower index “-I-” in 
w+ ). For Wr = w{x = d), one has the boundary condi¬ 
tion, d^ri;(a;)/da;^ = 0 at the ES {x = 0): 


e+[e{wr)\ -I- Wr(l + (iwr 


de+[e(u;)] 


dw 


= 0. (10) 


In Eq. (9), a « 0.5 —0.6 fm is the diffuseness parameter 

[7], 


a = 


C+pK+ 

mbl 


( 11 ) 


found from the asymptotic behavior of the particle den¬ 
sity, w ^ exp(—^/a) for large ^ a). 

As shown in Refs. [3, 5], the influence of the semi- 
classical h corrections (related to the ETF kinetic en¬ 
ergy) to w{x) is negligibly small everywhere, except for 
the quantum tail outside the nucleus (x^ 1). Therefore, 
all these corrections were neglected in Eq. (2). With 
a good convergence of the expansion of the e+[e(y)] in 
powers of 1 — y up to the leading quadratic term [3, 5] 
and small corrections in Eq. (4), e = (1 — y)^, one 
finds analytical solutions of Eq. (9) in terms of the al¬ 
gebraic, trigonometric and logarithmic functions [7]. For 
/3 = 0 (i.e. without SO terms), it simplifies to the solution 
w{x) = tanh^ [{x — xq)/2] for x < xq = 2arctanh(l/-\/3) 
and zero for x outside the nucleus {x > xq)- 

After simple transformations of the isovector Lagrange 
equation (Al) one similarly finds up to the leading term 
in a/H in the ES approximation for the isovector density, 
W-{x) = P-1 (jil), the equation and the boundary condi¬ 
tion (A3). The analytical solution w- = r(;cos[ 7 /’('u;)] can 
be obtained through the expansion (A5) of if) in powers 
of 


l{w) 


3e 

Csym 


Csym — ^ 



( 12 ) 


Expanding up to the second order in 7 one obtains (see 
Appendix A) 

W- = w cos [V'(i(')] « w ^1 — ^ + ...^ , (13) 

with 

/3Csym + 2 + Csy^L(l -b/3)/(3J) 

"" “ 2c,y™(l + /3) ’ ^ 

[see also the constant C 3 (Appendix A) at higher (third) 
order corrections]. The constant c [Eq. (15)] for the 
isovector solutions W-, Eq. (13), is modified with respect 
to Ref. [7] in two aspects. In addition to the L depen¬ 
dence there are also higher order terms from a nonlinear 
equation (A4) for tp{w) (Appendix A). Notice also that 
W- depends on L in second order in 7 but it is indepen¬ 
dent of K- at this order (Appendix A). 

In Fig. 1, the L dependence of the function w-{x) is 
shown within the total interval from L = 0 to L = 100 
MeV [20] and it is compared to that of the density w{x) 
for the SLy5* force as a typical example. As shown in 
Fig. 2 in a larger (logarithmic) scale, one observes no¬ 
table differences in the isovector densities W- derived 
from different Skyrme forces [37, 41] within the edge dif¬ 
fuseness. All these calculations have been done with the 
finite proper value of the slope parameter L. For SLy 
forces this value is taken from Ref. [44], for SGI I from 
Ref. [20] and for others from Ref. [41] (Table I). As shown 
below, this is in particular important for calculations of 
the neutron skin of nuclei. Notice that, with the precision 
of line thickness, our results are almost the same taking 
approximately L = 50 MeV for SLy5* and L = 60 MeV 
for SVsym32. Note also that, up to second order in the 
small parameter 7 , the isovector particle density W- in 
Eq. (13) does not depend on the symmetry-energy in¬ 
compressibility K-. The K- dependence appears only 
at higher (third) order terms in the expansion in 7 (Ap¬ 
pendix A). Therefore, as a first step of the iteration pro¬ 
cedure, it is possible to study first the main slope effects 
of L neglecting small P corrections to the isoscalar parti¬ 
cle density w+ (9) through e_|_ (4). Then, we may study 
more precisely the effect of the second derivatives K- 
taking into account higher order terms. 

We emphasize that the dimensionless densities, w{x) 
[see Eq. (9) and Ref. [7] ] and W-{x) (13), shown in 
Figs. 1 and 2, were obtained in leading ES approxima¬ 
tion {a/R <C 1) as functions of specific combinations of 
Skyrme force parameters like /3 and Csym [Eq. (12)] ac¬ 
counting for the L-dependence [Eq. (15)]. These densities 
are at the leading order in the leptodermous parameter 
a/R approximately universal functions, independent of 
the properties of the specihc nucleus. It yields largely 
the local density distributions in the normal-to-ES di¬ 
rection ^ with the correct asymptotic behavior outside 
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of the deformed ES layer at a/R <C 1, as is the case 
for semi-infinite nuclear matter. Therefore, at the domi¬ 
nating order, the particle densities w± are universal dis¬ 
tributions independent of the specific properties of nuclei 
while higher order corrections to the densities w± depend 
on the specific macroscopic properties of nuclei. 


For the energy surface coefficients (19), one ob¬ 
tains 

= 6C+pf7+/(roa), 

1 

J+= J dwy^w{l + I3w) e+[e(w)], (20) 

0 


IV. ISOVECTOR ENERGY AND STIFFNESS 


The nuclear energy E [Eq. (1)] in the improved ES ap¬ 
proximation (Appendix B) is split into volume and sur¬ 
face terms [7], 


E -by A +J{N - Z)^/A +Es . (16) 


For the surface energy Es one obtains 

As = 4+)+4") 


(17) 


with the isoscalar (-I-) and isovector (-) surface compo¬ 
nents: 


Ei^^ = 


S 

47rrg ’ 


(18) 


where S is the surface area of the ES, are the 
isoscalar (-I-) and isovector (—) surface-energy constants. 


4 ^ « 87rrgC± J dC (^1 + ^P+ 


dp± 


(19) 


These constants are proportional to the corresponding 
surface tension coefficients a± = /[A-ktI) through the 

solutions (9) and (13) for p±(^) which can be taken into 
account in leading order of a/R (Appendix B). These co¬ 
efficients cr± are the same as found in the expressions for 
the capillary pressures of the macroscopic boundary con¬ 
ditions [see Ref. [7] with new values e± modified by L and 
K- derivative corrections of Eqs. (4) and (7)]. Within 
the improved ES approximation where also higher order 
corrections in the small parameter a/R are taken into 
account, we derived in Ref. [7] equations for the nuclear 
surface itself (see also Refs. [2, 3, 5]). For more exact 
isoscalar and isovector particle densities we account for 
the main terms in the next order of the parameter a/R 
in the Lagrange equations (see Eq. (Al) for the isovec¬ 
tor and Refs. [2, 3, 5] for the isoscalar case). Multiplying 
these equations by dp-/d^ and integrating them over the 
ES in the normal-to-surface direction ^ and using the so¬ 
lutions for w±{x) up to the leading orders [Eqs. (9) and 
(13)], one arrives at the ES equations in the form of the 
macroscopic boundary conditions [2, 3, 5, 7, 30, 47-49]. 
They ensure equilibrium through the equivalence of the 
volume and surface (capillary) pressure variations. As 
shown in Ref. [7], the latter ones are proportional to the 
corresponding surface tension coefficients a±. 


=ks kg = 6p C„ J-/{roa), (21) 


J- = j dw 
0 

X ■[cos('0) + 


I w e+[e(w)] 
1 -I- jdw 


wsm{tp) 


Csym •\/l + fd 

i<‘V t 


[1 -I- 2c^{w)] 

2'y{w) 


(3w 


1 -b 


1 + /3 


'-sym 


6 ( 1 + / 3 ) 


-sym 


i(l + P) 

s 

1 

~ 2 


.( 22 ) 


For 7 and c, see Eqs. (12) and (15), respectively. Simple 
expressions for the constants in Eqs. (20) and (21) 
can be easily derived in terms of algebraic and trigono¬ 
metric functions by calculating explicitly integrals over w 
for the quadratic form of e+[e(w)] [Eqs. (B3) and (B5)]. 
Note that in these derivations, we neglected curvature 
terms and, being of the same order, shell corrections, 
which have been discarded from the very beginning. The 
isovector energy terms were obtained within the ES ap¬ 
proximation with high accuracy up to the product of two 
small quantities, P and (a/R)"^. 

According to the macroscopic theory [7-10], one may 
dehne the isovector stiffness Q with respect to the differ¬ 
ence Rn — Rp between the neutron and proton radii as a 
dimensionless collective variable r, 


= Qr^ 

T — (Rn Rp) /^O J 


Qt^S 

dTTPg ’ 


(23) 


where r is the relative neutron skin. Comparing this ex¬ 
pression to Eq. (18) for the isovector surface energy writ¬ 
ten through the isovector surface-energy constant ^ 
[Eq. (21)], one obtains 



(24) 


Defining the neutron and proton radii R„,p as positions 
of maxima of the neutron and proton density gradients, 
respectively, one obtains the neutron skin t (Ref. [7]), 


T = 


Sal 

’"OCgym 


g{wr), 


(25) 
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where 


= (l + /i)(3„. + l+4W 

+ 27 (1 + C 7 ) \cw - Csym (1 + 2c7)]} 


is taken at the ES value Wr [Eq. (10)]. Finally taking 
into account Eqs. (24) and (21), one arrives at 


^ _ MZ! - 

fcg’ 16 g‘^{wr)' 


(27) 


where JL and g{w) are given by Eqs. (22) and (26), re¬ 
spectively. Note that Q = —dJ'^/Akg was predicted in 
Refs. [8, 9] and therefore for v = 9/4 the first part of (27) 
which relates Q with the volume symmetry energy J and 
the isovector surface-energy constant kg, is identical to 
that used in Refs. [8-11, 19, 20]. However, in our deriva¬ 
tions V deviates from 9/4 and it is proportional to the 
function Jljg'^{wr)- This function depends significantly 
on the SO interaction parameter f3 but not too much on 
the specific Skyrme force (see Ref. [7] for details). 

Notice that the approximate universal functions w{x) 
[Eq. (9) and Ref. [7]] and W-{x) [Eq. (13)] can be used 
in the leading order of the ES approximation for calcu¬ 
lations of the surface energy coefficients [Eq. (19)] 
and the neutron skin r oc / [Eq. (25)]. As shown in 
Ref. [7] and in Appendix B, here only the particle den¬ 
sity distributions w{x) and w-{x) are needed within the 
surface layer through their derivatives [the lower limit of 
the integration over ^ in Eq. (19) can be approximately 
extended to —00 because there are no contributions from 
the internal volume region in the evaluation of the main 
surface terms of the pressure and energy]. Therefore, 
the surface symmetry-energy coefficient kg in Eqs. (21) 
and (B5) , the neutron skin r [Eq. (25)] and the isovector 
stiffness Q [Eq. (27)] can be approximated analytically in 
terms of functions of definite critical combinations of the 
Skyrme parameters like /3, Csym, cl, C- and parameters of 
infinite nuclear matter (6v, p, K+), also the symmetry 
energy constants J, L and K-. Thus, in the considered 
ES approximation, they do not depend on the specific 
properties of the nucleus (for instance, the neutron and 
proton numbers), the curvature and the deformation of 
the nuclear surface. 


V. THE FERMI-LIQUID DROPLET MODEL 


For IVDR calculations, the ELD model based on 
the linearized Landau-Vlasov equations for the isoscalar 
[(5/_^(r, p, t)] and isovector [6f_{r, p, t)] distribution func¬ 
tions can be used in phase space [30, 50, 51], 


dSU 

dt 



(28) 


Here e = p^/{2m’^) is the equilibrium quasiparticle en¬ 
ergy (p = IpI) and /{2m*^) is the Fermi en¬ 

ergy. The isotopic dependence of the Fermi momenta 
= Pp {1^ A) is given by a small parameter A = 
2{1 + Fq) 1/3 . The reason for having A is the differ¬ 
ence between the neutron and proton potential depths 
from the Coulomb interaction. The isotropic isoscalar 
Fq and isovector F/ Landau interaction constants are re¬ 
lated to the isoscalar in-compressibility K = 06^(1 -I- Fq) 
and the volume symmetry energy J = 2ej,,(l -I- F/)/3 
constants of nuclear matter, respectively. The effective 
masses = to(1-|-Ei/ 3) and mZ_ = to( 1-|-E(/3) are de¬ 
termined in terms of the nucleon mass m by anisotropic 
Landau constants Fi and F[. Equations (28) are cou¬ 
pled by the dynamical variation of the quasiparticles’ 
selfconsistent interaction 5V± with respect to the equilib¬ 
rium value p^ l(2m*p). The time-dependent external field 
V/xt oc exp(—zwt) is periodic with a frequency w. For sim¬ 
plicity, the collision term 6St± is calculated within the 
relaxation time T(w) approximation accounting for the 
retardation effects from the energy-dependent self-energy 
beyond the mean field approach, T = 47r^7o/(fia;)^ with 
the parameter To oc [see Eq. (80) of Ref. [51] at 

zero temperature and also Ref. [30]]. 

The solutions of Eq. (28) are related to the dynamic 
multipole particle-density variations, Sp^{r,t) oc Y\o{f), 
where lAo(r) are the spherical harmonics and f = r/r. 
These solutions can be found in terms of the superposi¬ 
tion of plane waves over the angle of a wave vector q , 

Sf±iP,r,t) = J dQqVAo(T) Sf±{p,q,uj) 

X exp [—i {cot — qr)], (29) 

where i5/±(p, q, w) is the Fourier transform of the dis¬ 
tribution function. The time-dependence (29) is peri¬ 
odic as the external field is also periodic with the 
same frequency w = p^s^q/m’^ where s’*" = s, and 

s~ = s (fVZ/A^)The factor accounts 

for conserving the position of the mass center for the 
isovector vibrations [53]. The sound velocity s can be 
found from the dispersion equations [30]. The two solu¬ 
tions Sn with n = 1,2 are functions of the Landau interac¬ 
tion constants and ujF- Due to the symmetry interaction 
coupling the “out-of-phase” particle-density vibrations of 
the Si mode involve the “in-phase” mode S2 inside of the 
nucleus. 

For small isovector and isoscalar multipole ES-radius 
vibrations of the finite neutron and proton Fermi- 
liquid drops around the spherical nuclear shape, one 
has SR±{t) = Rag{t)Y\o{f) with a small time- 
dependent amplitudes a^{t) = a^exp{—iujt). The 
macroscopic boundary conditions (surface continuity and 
force-equilibrium equations) at the ES are given by Refs. 
[7, 30, 51]: 

uf = RcegYxoif), 

r=R 

= a^Pg Yxoif). ( 30 ) 

r—R 
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The left hand sides of these equations are the radial com¬ 
ponents of the mean-velocity field u = j/p (j is the cur¬ 
rent density) and the momentum flux tensor (511^^ defined 
both through the moments of Sf{r,p,t) in momentum 
space [30, 51]. The right-hand sides of Eq. (30) are the 
ES velocities and capillary pressures. These pressures are 
proportional to the isoscalar and isovector surface-energy 
constants 6^ in Eq. (19), 

Pt = (31) 

where V+ = (A — 1)(A + 2)/2 , V- = 1 . The coeffi¬ 
cients are essentially determined by the constants C± 
[Eq. (8)] of the energy density (2) in front of its gradient 
density terms. The conservation of the center of mass 
is taken into account in the derivations of the second 
boundary conditions (30) [30, 51]. Therefore, one has a 
dynamical equilibrium of the forces acting at the ES. 


VI. TRANSITION DENSITY AND NUCLEAR 
RESPONSE 


The response function, x±(w), is defined as a linear 
reaction to the external single particle field F(r) with 
the frequency cu. For convenience, we may consider this 
field in terms of a similar superposition of plane waves 
(29) as Sf± [30, 51]. In the following, we will consider the 
long wave-length limit with 

^ q;^’“ is the amplitude 

and w is the frequency of the external field {rjo = -|-0). 
In this limit, the one-body operator F{r) becomes the 
standard multipole operator, F{r) = r^Y\o{f) for A > 1. 
The response function x±(w) is expressed through the 
Fourier transform of the transition density (r) as 

X±(w) = - y dr ^(r) p±(r)/a^^’“. ( 32 ) 


The transition density p^ (r) is obtained through the dy¬ 
namical part of the particle density dp±{r, t) in a macro¬ 
scopic model in terms of solutions Sf±{r,p,t) of the 
Landau-Vlasov equations (28) with the boundary con¬ 
ditions (30) as the same superpositions of plane waves 
(29) [30]: Sp-{r,t) = p asP-{x) Tio(r) where 


P-{^) 


qR 

i'MR) 


(k) w{x) + 


9v dw- ' 

9s dx J ’ 


(33) 


9v = [ dw 

J l-w 

0 

Wq 

9s = j dw [1 + 0(7 ^(w))] , 


0 

K = tin 


U/ , - 

1 + ;^2 :(w) 


(34) 

(35) 

(36) 


Ko = qR- The first term in (33), proportional to the 
dimensionless isoscalar density w(x) (in units of p) ac¬ 
counts for volume density vibrations [Eq. (9)]. The sec¬ 
ond term oc dw-/dx, where W- is a dimensionless isovec¬ 
tor density (in units of p/) corresponds to the density 
variations from a shift of the ES [Eq. (13)]. The particle 
number and the center-of-mass position are conserved, 
and j\{K) and j'xiK) are the spherical Bessel functions 
and their derivatives. The upper integration limit iq, in 
Eqs. (34) and (35) is defined as the root of a transcendent 
equation x{wo) + R/a = 0. As shown in Appendix A, the 
SO and L dependent density W-{x) is of the same order 
as w{x). The dependencies of W-{x) on different Skyrme 
force parameters, mostly the isovector gradient-term con¬ 
stant C-, the SO parameter /3, and the derivative of the 
volume symmetry energy L are the main reasons for the 
different values of the neutron skin. 

With the help of the boundary conditions (30), one can 
derive the response function (32) [30], 

(uj - 10 , 

n n ^ ' 

(37) 

with uj = ppSnKo /{m*R) {m*_ « m\ = 

m*). This response function describes two modes, the 
main (n = 1) IVDR and its satellite (n = 2) as related 
to the out-of-phase si and in-phase S2 sound velocities 
which are excited in the nuclear volume, respectively. 
We assume here that the “main” peak exhausts mostly 
the energy weighted sum rule (EWSR) and the “satel¬ 
lite” corresponds to a much smaller part of the EWSR as 
proportional to the asymmetry parameter, / <C 1. This 
two-peak structure is from the coupling of the isovector 
and isoscalar density-volume vibrations because of the 
neutron and proton quasiparticle interaction 5V± in Eq. 
(28). Therefore, one takes into account an admixture of 
the isoscalar mode to the isovector IVDR excitation. The 
wave numbers q = Ko/R of the lowest poles (n = I, 2) in 
the response function (37) are determined by the secular 
equation. 


= j'xino)- 


“iCpKoC^l^ 
2&sAi/3 1 


Ja(«;o) +4"^j"(Ko) 


= 0 . 


(38) 

The width of an IVDR peak T in (37) corresponds to 
an imaginary part of the pole having its origin in the 
collision term 5St± of the Landau-Vlasov equation. At 
this pole, for the relaxation time one has 


Tn = ATT^%/{hWnY (39) 

with an A-dependent constant, To oc A~^/^. For the 
amplitudes one has oc A”“^. The complete expres¬ 
sions for the amplitudes and the constants are 
given in Refs. [30, 51]. Assuming a small value of A, 
one may call the n = 2 mode a “satellite” to the “main” 
n = 1 peak. On the other hand, other factors such as a 
collisional relaxation time, the surface symmetry-energy 



7 


constant bg, and the particle number A lead sometimes 
to a re-distribution of the EWSR values among these 
two IVDR peaks. The slope L dependence of the tran¬ 
sition densities p‘^{x) [Eq. (33)] and the strength of the 
response function, 

S'(w) = ImxA(w)/7r (40) 

[Eq. (37)] has its origin in the symmetry-energy coeffi¬ 
cient [Eqs. (21), (22), (15), (4) and (6)]. Thus, one 
may evaluate the EWSR sum rule contribution of the 
nth peak by integration over the region around the peak 
energy = hiOn, 

J doj uj Sn{(jj). (41) 

In accordance with the time-dependent HE approaches 
based on the Skyrme forces, (see, for instance, [24, 25, 
28]), we may expect that the energies of the satellite 
resonances in the IVDR and ISDR channels can be 
close. Therefore, we may calculate separately the neu¬ 
tron, Pn{x), and proton, Pp{x), transition densities for 
the satellite by calculating the isovector and isoscalar 
transition densities at the same energy E 2 and in the 
same units as p±, 

, p%{x)+p‘t{x) ^ p+{x)-p‘t{x) 

Pn{x) = — -^-, Pp{x) = - - -. 

(42) 

VII. DISCUSSION OF THE RESULTS 

In Table II we show the isovector surface energy coef¬ 
ficient kg [Eq. (21)], the stiffness parameter Q [Eq. (27)], 
its constant v and the neutron skin r [Eq. (25)] for many 
more Skyrme forces than discussed in [36]. They are ob¬ 
tained within the ES approximation with the quadratic 
expansion for e+[e(t(;)] and neglecting the P slope cor¬ 
rections, for several Skyrme forces [37, 41] whose pa¬ 
rameters are presented in Table I. Also shown are the 
quantities kso, 1 ^ 0 , Qo and tq neglecting the slope correc¬ 
tions {L = 0,Ar_ = 0). This is in addition to results of 
Ref. [7] where another important dependence on the SO 
interaction measured by /3 was presented. In contrast 
to a fairly good agreement for the analytical isoscalar 
surface-energy constant h^g'' (20) as shown in Ref. [7] 
and references cited therein, the isovector energy coef¬ 
ficient kg is more sensitive to the choice of the Skyrme 
forces than the isoscalar one b^g'^ [Eq. (20) and Ref. [5]]. 
The modulus of kg is significantly larger for most of the 
Skyrme forces SLy... [37] and SV... [41] than for the other 
ones. However, the L dependence of kg is somewhat 
small in these forces (cf. the first two rows of Table H) 
as it should be for a small parameter e of the symmetry- 
energy density expansion (6). For SLy and SV forces, the 
stiffnesses Q are correspondingly significantly smaller in 
absolute value being closer to the well-known empirical 


values Q ~ 30 — 35 MeV [9-11] obtained by Swiatecki 
and collaborators. Note that the isovector stiffness Q is 
even much more sensitive to the parametrization of the 
Skyrme force and to the slope parameter L than the con¬ 
stants kg. In Ref. [7], we studied the hydrodynamical 
results for Q as compared to the ELD model for the av¬ 
eraged properties of the giant IVDR (IVGDR) at zero 
slope L = 0. The IVDR structure in terms of the two 
(main and satellite) peaks was discussed earlier in Refs. 
[31] at L = 0 in some magic nuclei with a large neu¬ 
tron excess within the semiclassical ELD model based 
on the effective surface approach. For the comparison 
with experimental data and other theoretical results we 
present in Table H (row 9 and 11) a small L dependence of 
the IVGDR energy parameter D = EivgdrA^/^, where 
Aivgdr = [EiSi + E 2 S 2 ] / [S'! -k S'a] is the IVGDR en¬ 
ergy for the isotope ^^^Sn [ S'„ = S'(w„) , see also Eq. (40) 
for the definition of the strength 5'(u;)]. A more precise 
reproduction of the A-dependence of the IVGDR energy 
parameter D for finite values of L (see the last three 
rows for several isotopes) might determine more consis¬ 
tent values of Q, but, at present, it seems to be beyond 
the accuracy of both the hydrodynamical and the ELD 
models. The IVGDR energies obtained by the semiclas¬ 
sical Landau-Vlasov equation (28) with the macroscopic 
boundary conditions (30) of the ELD model (Ref. [7]) are 
also basically insensitive to the isovector surface energy 
constant ks [6, 7, 31, 36]. They are in a good agreement 
with the experimental data, and do not depend much on 
the Skyrme forces even if we take into account the slope 
symmetry-energy parameter L (last three rows in Table 
H). 

More realistic self-consistent HE calculations taking 
into account the Coulomb interaction, the surface cur¬ 
vature, and quantum shell effects have led to larger val¬ 
ues of (5 ~ 30 — 80 MeV [4, 20]. For larger Q (see Ta¬ 
ble H) the fundamental parameter (9 J/4(5)A“^/^ of the 
LDM expansion in Ref. [8] is really small for A> 40, and 
therefore, the results obtained using the leptodermous 
expansion are better justified. 

An investigation within the approach presented in Sec. 
V shows that the IVDR strength is split into a main 
peak which exhausts an essential part of the EWSR in¬ 
dependent of the model and a satellite peak with a much 
smaller contribution to this quantity. Focusing on a much 
more sensitive ks dependence of the IVDR satellite res¬ 
onances, one may take now into account the slope L de¬ 
pendence of the symmetry-energy density per particle 
(6) (Refs. [22-25] and [31]). The total IVDR strength 
function being the sum of the “out-of-phase” n = I and 
“in-phase” n = 2 modes for the isovector- and isoscalar¬ 
like volume particle density vibrations, respectively (solid 
lines in Figs. 3 and 4 for the zero L and dotted and dashed 
lines for the finite L) has a somewhat remarkable shape 
asymmetry [31, 36]. For SLy5* (Fig. 3) and for SVsym32 
(Fig. 4) one has the “in-phase” satellite to the right of the 
main “out-of-phase” peak, cf. with the traditional PDR 
to the left of the main one. An enhancement to the left 



of the main peak for SLy5* is from increasing the “out- 
of-phase” strength (red solid and magenta rare dotted 
curves, Fig. 3) at small energies because of the appear¬ 
ance of a peak at the energy about a few MeV, in contrast 
to the SVsym32 case. The semiclassical FLD model cal¬ 
culations at the lowest h order should be improved here, 
for instance by taking into account the quantum effects 
like shell corrections within more general periodic-orbit 
theory [51, 52]. In the nucleus ^^^Sn the IVDR energies of 
the two peaks do not change much with L in both cases: 
El = 17 MeV, E 2 = 20 MeV for SLy5* (Fig. 3) and 
El = 15 MeV, E 2 = 18 MeV for SVsym32 (Fig. 4). We 
find only an essential re-distribution of the EWSR con¬ 
tributions (normalized to 100% for the EWSR sum of the 
main and satellite peaks) [Eq. (41) for This is from 
a significant enhancement of the main “out-of-phase” 
peak with increasing L, s[^'^ = 89% and 52^^ = 11% for 
SLy5* (Eig. 3) and more pronounced EWSR distribution 
S'W = 76% and = 24% for SVsym32 (Eig. 4) [cf. 
with the corresponding L = 0 results: s[^'^ = 88% and 
= 12% for SLy5* and ^ = 73% and = 27% 
for SVsym32]. 

Eigures 5 and 6 show more systematic study for several 
isotopes and for the chain of the Sn isotopes, respectively. 
In Fig. 6, we compare the results of our calculations 
with the experimental data. The latter were obtained 
by fitting the experimental strength curve for a given al¬ 
most spherical Sn isotope by the two Lorenzian oscillator 
strength functions as described in Refs. [30, 51]. It is al¬ 
ways possible in the case of the asymmetric shapes of the 
strength curves with usual enhancement on the right of 
the main peak, even in the case when the satellite cannot 
be distiguished transparently well from the main peak in 
almost spherical nuclei (unlike the clear shoulders for the 
IVDRs in deformed ones). Each of these functions has 
three fitting parameters such as the inertia, stiffness and 
width of the peak. We found somewhat good agreement 
of our ETE ES results with these experimental data for 
the energies, ratio of the strengths at the satellite to the 
main modes, and the EWSR contributions. 

More precise L-dependent calculations change essen¬ 
tially the IVDR strength distribution for the SV forces 
because of the smaller Cgym value as compared to other 
Skyrme interactions (see Table I). For ^°®Pb one obtains 
El = 15 MeV, = 91% for the main peak and E 2 = 17 
MeV, S'2^^ = 9% the satellite for SLy5*; and Ei = 13 
MeV, = 83% for the main peak and E 2 = 16 MeV, 
S'T) = 17% tiie satellite for SVsym32 forces. These cal¬ 
culations are qualitatively in agreement with the experi¬ 
mental results: Ei = 13 MeV, = 98% for the main 
peak and E 2 = 17 MeV, = 2% the satellite. Dis¬ 
crepancies might be related to the strong shell effects in 
this stable double magic nucleus which are neglected in 
the ETF ES approach. 

Decreasing the relaxation time T by a factor of about 
1.5 almost does not change the IVDR strength structure. 


However, we found a strong dependence on the relaxation 
time T in a wider region of E values. The “in-phase” 
strength component with a wide maximum does not de¬ 
pend much on the Skyrme force [37, 41, 43], the slope 
parameter L, and the relaxation time 7”. We found also 
a regular change of the IVDR strength for different dou¬ 
ble magic isotopes (Eig. 5). In addition to a big change 
for the energy (mainly because of Ei) and the strength 
[<S'i(u;)], one also obtains more asymmetry for ®®Ni than 
for the other isotopes. Calculations for nuclei with differ¬ 
ent mass A were performed with the relaxation time E 
[Eq. (39)] where To = 7opb(208/H)^/^ with the parame¬ 
ter ToPb = 300 MeV^-s derived from the IVGDR width of 
^°®Pb, in agreement with experimental data for the av¬ 
eraged A dependence of the IVGDR widths (oc 
In this way the IVDR relaxation time En becomes larger 
with increasing A as and at the same time, the 

height of peaks decreases. The L corrections are also 
changing much in the same scale of all three nuclei. 

The essential parameter of the Skyrme HF approach 
leading to the significant differences in the ks and Q val¬ 
ues is the constant C_ [Eq. (2) and Table I]. Indeed, C_ 
is the key quantity in the expression for Q [Eq. (27)] 
and the isovector surface-energy constant kg [or bg ^, Eq. 
(21)], because Q oc 1/ks oc 1/C_ and kg oc C_ [7]. Con¬ 
cerning kg and the IVDR strength structure, this is even 
more important than the L dependence although the lat¬ 
ter changes significantly the isovector stiffness Q and the 
neutron skin r. As seen in Table I, the constant C_ is 
very different in absolute value and in sign for differ¬ 
ent Skyrme forces whereas C+ is almost constant (Table 
I). The isoscalar energy density constant 6^"^ is propor¬ 
tional to C+ [Eq. (20)], in contrast to the isovector one. 
All Skyrme parameters are fitted to the well-known ex¬ 
perimental value 6^"^ = 17—19 MeV while there are 
so far no clear experiments which would determine kg 
well enough because the mean energies of the IVGDR 
(main peaks) do not depend very much on kg for differ¬ 
ent Skyrme forces (the last three rows of Table II). Per¬ 
haps, the low-lying isovector collective states are more 
sensitive but at the present time there is no careful sys¬ 
tematic study of their kg dependence. Another reason 
for so different ks and Q values might be because of dif¬ 
ficulties in deducing kg directly from the HF calculations 
because of the curvature and quantum effects. In this re¬ 
spect, the semi-infinite Fermi system with a hard plane 
wall might be more adequate for the comparison of the 
HF theory and the ETF effective surface approach. We 
have also to go far away from the nuclear stability line 
to subtract uniquely the coefficient kg in the dependence 
of (X E = {N — Z)‘^IA‘^, according to Eq. (21). Eor 
exotic nuclei one has more problems to derive kg from 
the experimental data with enough precision. Note that, 
for studying the IVDR structure, the quantity kg is more 
fundamental than the isovector stiffness Q because of the 
direct relation to the tension coefficient (t_ of the isovec¬ 
tor capillary pressure. Therefore, it is simpler to analyze 
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the experimental data for the IVGDR within the macro¬ 
scopic HD or FLD models in terms of the constant kg. 
The quantity Q involves also the ES approximation for 
the description of the nuclear edge through the neutron 
skin T in Eq. (24). The L dependence of the neutron 
skin r is essential but not so dramatic in the case of SLy 
and SV forces (Table II), except for the SVmasOS forces 
with the effective mass 0.8. The precision of such a de¬ 
scription depends more on the specific nuclear models 
[19, 20, 27]. On the other hand, the neutron skin thick¬ 
ness r, like the stiffness Q, is interesting in many aspects 
for an investigation of exotic nuclei, in particular, in nu¬ 
clear astrophysics. 

We emphasize that for specific Skyrme forces there ex¬ 
ists an abnormal behavior of the isovector surface con¬ 
stants kg and Q. It is related to the fundamental con¬ 
stant C- of the energy density (2) but not to the deriva¬ 
tive symmetry-energy density corrections. For the pa¬ 
rameter set T6 (C_ = 0) one finds kg = 0 (Ref. [7]). 
Therefore, according to Eq. (27), the value of Q diverges 
{ly is almost independent from C_ for SLy and SV forces; 
Table II and Refs. [7, 31, 36]). The isovector gradient 
terms which are important for the consistent derivations 
within the ES approach are also not included (C_ = 0) 
into the symmetry energy density in Refs. [15, 17, 18]. In 
relativistic investigations [12, 13] of the structure of the 
IVGR distributions, the dependence of these quantities 
on the derivative terms has not been investigated so far. 
It therefore remains an interesting task for the future to 
apply similar semiclassical methods such as the ES ap¬ 
proximation used here also in relativistic models. More¬ 
over, for RATP [37] and SV [41] (like for Ski) Skyrme 
forces, the isovector stiffness Q is even negative as C_ >0 
{kg > 0) in contrast to other Skyrme forces. This would 
lead to an instability of the vibration of the neutron skin. 

Table II shows also the coefficients ly of Eq. (27) for 
the isovector stiffness Q. They are almost constant for 
all SLy and SV Skyrme forces, in contrast to other forces 
[7]. However, these constants ly, being sensitive to the SO 
(/3) dependence through Eqs. (26), (25) and (22), change 
also with L (Table H). As compared to 9/4 suggested in 
Ref. [8], they are significantly smaller in magnitude for 
most of the Skyrme forces. 

In Fig. 7 we show, in the case of the Skyrme forces 
SLy5* and SVsym32, the transition densities p‘^{x) of 
Eq. (33) for the “out-”of-phase (-) and the “in-” phase 
(-I-) modes of the volume vibrations at the excitation en¬ 
ergy E 2 of the satellite. The transition densities were 
not presented in our preliminary publication [36]. These 
are the key quantities for the calculation of the IVDR 
strengths, according to Eq. (32). The L dependence 
is somewhat small, slightly notable mostly near the ES 
(ja;j 1). Erom Eig. 8, one finds a remarkable neutron vs 
proton excess near the nuclear edge for the same forces, 
which is however, very slightly depending on the slope 
parameter L. A small dependence of the transition den¬ 
sities on L comes through the symmetry-energy constant 
kg which is almost the same in modulus for these forces. 


We did not find a dramatic change of the transition den¬ 
sities with the sign of kg. Therefore, there is a weak 
sensitivity of the transition densities on L through the 
energy E 2 . We would have expected a stronger influence 
of the sign of kg on the vibrations of the neutron skin 
rather than on the IVDR. This different sign leads to the 
opposite, stable and unstable, neutron skin vibrations. 
One observes also other differences between the upper 
(SLy5*) and the lower (SVsym32) panels in both figures: 
We find a redistribution of the surface-to-volume contri¬ 
butions of the transition densities for these two modes. 
Again, as in Eigs. 9 and 10, one finds a considerable 
change of the neutron-proton transition densities for the 
different isotopes for SLy5* and SVsym32. 

The last figure shows theoretical (Fig. II) evaluations 
of the neutron skin. Figure II shows the absolute val¬ 
ues of the skin obtained from t// multiplying the mean- 
square evaluations of the nuclear radii by the factor a/3/5 
for an easy comparison with experimental data given in 
[27]. For ^°®Pb, one finds that the experimental val¬ 
ues = 0.12 - 0.14 fm in Ref. [27] (0.156;[!/2i 

fm [59]) are in good agreement with our calculations 
~ 0-10 ~ 0-13 fm within the ES approximation 
(the limits show values from SLy5* to SVsym32). For the 
isotope ^^"^Sn one obtains Ar^p®"” « 0.09 — 0.12 fm, also 
in good agreement with experimental results. For the 
isotope ^^^Sn, we predict the value Ar^jp®®” « 0.11 — 0.15 
fm. Similarly, for ®°Ni and ®®Ni, one finds « 

0.03 — 0.04 (like in Ref. [27]) and 0.08 — 0.11 fm, respec¬ 
tively. 

VIII. CONCLUSIONS 

The slope parameter L was taken into account in the 
leading ES approximation to derive simple analytical ex¬ 
pressions for the isovector particle densities and energies. 
These expressions were used for calculations of the sur¬ 
face symmetry energy, the neutron skin thickness, and 
the isovector stiffness coefficients as functions of L. For 
the derivation of the surface symmetry energy and its de¬ 
pendence on the particle density we have to include main 
higher order terms in the parameter a/R. These terms 
depend on the well-known parameters of the Skyrme 
forces. Results for the isovector surface-energy constant 
fcs, the neutron skin thickness r and the stiffness Q de¬ 
pend in a sensitive way on the parameters of the Skyrme 
functional (especially on the parameter C_) in the gradi¬ 
ent terms of the density in the surface symmetry energy 
[see Eq. (2)]. The isovector constants fcs, r and Q depend 
also essentially on the slope parameter L, in addition to 
the SO interaction constant /3. Eor all Skyrme forces, 
the isovector stiffness constants Q are significantly larger 
than those obtained in earlier investigations. However, 
taking into account their L-dependence they come closer 
to the empirical data. It influences more on the isovector 
stiffness Q and on the neutron skin t, than on the surface 
symmetry-energy constant kg. The mean IVGDR ener- 
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gies and sum rules calculated in the macroscopic models 
like the FLD model [6, 30] in Table II are in fairly good 
agreement with the experimental data for most of the 
ks values. As compared with the experimental data and 
other recent theoretical works, we found a somewhat rea¬ 
sonable two-peak structure of the IVDR strength within 
the FLD model. According to our results for the neutron 
and proton transition densities [Figs. 8-10], we may inter¬ 
pret semiclassically the IVDR satellites as some kind of 
pygmy resonances, in addition to the traditional studies 
[12-14, 22-25, 32-35]. Their energies, sum rules and n-p 
transition densities obtained analytically within the semi- 
classical FLD approximation are sensitive to the surface 
symmetry-energy constant and the slope parameter L. 
Therefore, their comparison with the experimental data 
can be used for the evaluation of fcg and L. It seems 
helpful to describe them in terms of only few critical pa¬ 
rameters, like kg and L. 

For further perspectives, it would be worthwhile to 
apply our results to calculations of the IVDR strength 
structure within the FLD model [30] in a more system¬ 
atic way. In this respect it is also interesting that the 
low-lying collective isovector states are expected to be 
even more sensitive to the values of ks within the peri¬ 
odic orbit theory [52, 60, 61]. More general problems 
of classical and quantum chaos in terms of the level 
statistics and Poincare and Lyapunov exponents (Ref. 
[62] and references therein) might lead to progress in 
studying the fundamental properties of collective dynam¬ 
ics like nuclear fission within the Swiatecki-Strutinsky 
macroscopic-microscopic model. Our approach is helpful 
also for further study of the effects in the surface sym¬ 
metry energy because it gives analytical universal expres¬ 
sions for the constants kg, r and Q as functions of the 
slope parameter L which do not depend on specific prop¬ 
erties of nuclei as they are directly connected with a few 
critical parameters of the Skyrme interaction without any 
fitting. 
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The Lagrange equation for the variations of the isovec¬ 
tor particle density p_ is given in the local coordinates 
C,.by[5,7] 


2C_ 


d^P- 

de 


2C-n 




- — [p+£- {p+,p-)] + A- = 0, 
dp_ 


(Al) 


where H is the mean curvature of the ES, A_ is the ES 
correction to the isovector chemical potential. Up to the 
leading terms in a small parameter a/R one gets from 
Eq. (AI) 


2C_ 




— [p+£-{p+,P-)] = 0 • 
dp- 


(A2) 


We neglected here the higher order terms proportional to 
the first derivatives of the particle density p_ with respect 
to ^ and the surface correction to the isovector chemical 
potential in Eq. (AI) (Refs. [2, 3] for the isoscalar case). 
For the dimensionless isovector density W- = p-/{j)I) 
one finds after simple transformations the following equa¬ 
tion and the boundary condition in the form 


dW- _ / 5sym(£)(l + Pw) I ^ wt 

dw e[e{w)] V 

W- (w = 1) = I , 


(A3) 


where P is the SO parameter defined below Eq. (9), 
tSsym ~ ^sym/Ji Csym Is defined in Eq. (12) and tSsyni(£) 
in Eq. (6). The above equation determines the isovector 
density W- as a function of the isoscalar one w{x) [Eq. 
(9)]. In the quadratic approximation for e[e(t(;)] [up to 
a small asymmetry correction proportional to in Eq. 
(4)], one finds an explicit analytical expression in terms 
of elementary functions [7]. Substituting w- = w cosf/' 
into Eq. (A3), and taking the approximation e = (1—rc)^, 
one has the following first order differential equation for 
a new function ipiw)'. 
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Appendix A: SOLUTIONS OF THE 
ISOVECTOR LAGRANGE EQUATION 


t(;(I — w) . , dib 

-sm^ -— = 

Csym d'W 

- COS'(/>, ' 4 !{w 

Csym 


\JSsym{(){l + Pw) simp 
= I) = 0. (A4) 


The boundary condition for this equation is related 
to that of Eq. (A3) for w-{w). This equation looks 
more complicated because of the trigonometric nonlinear 
terms. However, it allows one to obtain simple approx¬ 
imate analytical solutions within standard perturbation 
theory. Indeed, according to Eqs. (A3) and (9), where 
we do not have an explicit x-dependence, we note that 
W- is mainly a sharply decreasing function of x through 
w{x) within a small diffuseness region of the order of one 
in dimensionless units (Figs. I and 2). Thus, we may 
find approximate solutions to Eq. (A4) with its bound¬ 
ary condition in terms of a power expansion of a new 
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function ipij) in terms of a new small argument 7 [Eq. coordinates r in the local coordinate system 7 (see Fig. 
(12)] 1) in the ES approximation, one finds 


^( 7 ) = = '^Cn 7”(w) , (A5) 

n—0 


where the coefficients c„ and 7 are defined in Eq. (12). 
Substituting the power series (A5) into Eq. (A4), one 
expands first the trigonometric functions into a power 
series of 7 according to the boundary condition in Eq. 
(A4). As usual, using standard perturbation theory, we 
obtain a system of algebraic equations for the coefficients 
c„ [Eq. (A5)] by equating coefficients from both sides of 
Eq. (A4) with the same powers of 7 . This simple proce¬ 
dure leads to a system of algebraic recurrence relations 
which determine the coefficients c„ as functions of the 
parameters (3 and Csym of Eq. (A4), 


Co = 0 , Cl = 

Cl 


1 


C 2 = 


2Csym(l + P) 


C 3 = -Cl 

1 


c?-3 


yiT?’ 

Pc: 

C1C2 C 2 Cgyni 


.2 

'-sym 


+ 2 + ^Cgyjjj(l -f /3) 


‘-"sym 


2 ci 


Pcf 


L 


- 


K 

4 , —'^sym 

^ 36J 


(A 6 ) 


12 J V 1 6j) 


and so on. In particular, up to second order in 7 , we 
derive analytical solutions as functions of /3, Cgym, J and 
L in an explicitly closed form: 

^(7) = 7 (ci + C27), Cl = (A 7 ) 

/^Cgym + 2 + Acsy^(l-I-/3)/(3J) 

Thus, using the standard perturbation expansion method 
of solving tpij) in terms of the power series of the 7 up 
to 7 ^, one obtains the quadratic expansion of 'ip{w) [Eq. 
(13)] with c = C 2 /C 1 . Notice that one finds a good con¬ 
vergence of the power expansion of 'ip{j(w)) (A7) in ^{w) 
for w-{x) at the second order in 'y(w) because of values 
of Csym larger one for all Skyrme forces presented in Table 
I [Eq. (12) for Cgym]- 


Appendix B: Derivations of the surface energy 
and its coefficients 


For the calculation of the surface energy components 
of the energy E in Eq. (1) within the same improved 
ES approximation as described above in Appendix A we 
first separate the volume terms related to the first two 
terms of Eq. (2) for the energy density £ per particle. 
Other terms of the energy density p£{p+,p-) in Eq. (2) 
lead to the surface components Eg [Eq. (18)], as they are 
concentrated near the ES. Integrating the energy density 
p£ per unit of the volume [see Eq. (2)] over the spatial 


= <pdS / dC 


C±(Vp±) + p+e±{p+,p-) 


cr± S, 


(Bl) 


where — a (Refs. [2, 3, 5]). The local coordi¬ 

nates rj were used because the integral over ^ converges 
rapidly within the ES layer which is effectively taken for 
1^1 a. Therefore again, we may extend formally 
to —00 in the first (internal) integral taken over the ES 
in the normal direction ^ in Eq. (BI). Then, the sec¬ 
ond integration is performed over the closed surface of 
the ES. The integrand over ^ contains terms of the or¬ 
der of (p/a)^ oc {R/aY like the ones of the leading order 
in the first equation of Ref. [7]. However, the integra¬ 
tion is effectively performed over the edge region of the 
order of a that leads to the additional smallness propor¬ 
tional to a/R like in Appendix A. At this leading order 
the 7 dependence of the internal integrand can be ne¬ 
glected. Moreover, from the Lagrange equations [see Eq. 
(A2) for the isovector case] at this order one can realize 
that terms without the particle density gradients in Eq. 
(Bl) are equivalent to the gradient terms. Therefore, 
for the calculation of the internal integral we may ap¬ 
proximately reduce the integrand over ^ to derivatives of 
the universal particle densities of the leading order p±{^) 
in ^ using C± {Vp±f + p+e± (p+,p-) « 2C±{dp±/d^f 
[see Eqs. (9) and (13) for w±(a;)] . We emphasize that 
the isovector gradient terms are obviously important for 
these calculations. Taking the integral over ^ within the 
infinite integration region (—00 < ^ < 00 ) out of the in¬ 
tegral over the ES (dS') we are left with the integral over 
the ES itself that is the surface area S. Thus, we arrive 
finally at the right-hand side of Eq. (Bl) with the sur¬ 
face tension coefficient a± = /(AirrQ) [ see Eq. (19) 

for 6 g^^]. 

Using now the quadratic approximation e[e{w)] = (1 — 
wY in Eq. (19) for iV- = 0) one obtains (for /3 < 0, 
see Table I) 


=&-pC± ^ 

roa 


(B2) 


where 


1 

J+ = j dw a/w(1 -I- Pw) (1 - w) 


(B3) 


24(-/3)5/2 


+ P) + J+'’ arcsin a/^ 


with 


jyW + + (B4) 
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For the isovector energy constant J- one finds with 

= 105-4/3{95 + 75C + /3[119 + 10C(19 + 6C) 

+ 8/32 (1 + 10C(1 + 0) + 8C (5C(3 + 20 - 6)] } , 

J-\0 = 15 {7+ 2/3 [5(3 + 20+8/3(1 + 0 

X (3 + C + 2/3(1 + 0)]}. (B6) 


J- 


1 

^ ^ p j dw a/ w(l + /3w)(l - w)(l + c"f{w)Y 
0 




1920(1 +/3)(-/3)9/2 1 
+ 


(^) 0-/3(l + /3) 

^ Cgym ^ \/~/^ J (B5) 


These equations determine explicitly the analytical ex¬ 
pressions for the isoscalar and isovector ( 6 ^ en¬ 

ergy constants in terms of the Skyrme force parameters; 
see Eq. (15) for c and Eq. (12) for Csym and 7 (w). For the 
limit /3 —>■ 0 one has from Eqs. (B3) and (B5) J± 4/15. 
With Eqs. (25) and (26) one arrives also at the explicit 
analytical expression for the isovector stiffness Q as a 
function of C_ and /3. In the limit C_ —>■ 0 one obtains 
kg ^ 0 and Q —>■ 00 because of the finite limit of the 
argument Csym/c —>■ 2(1 + /3)/[/3 + (1 + /3)T/(3J)] of the 
function J7_ in Eq. (B5) [see also Eqs. (13) for c and Eq. 
(12) for Csym]- 
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SkM* 

SGII 

SLy5 

SLy5* 

SLy6 

SLy7 

SVsym28 

SVsym32 

SVmas08 

SVK226 

SVkap02 

P (fm-^) 

0.16 

0.16 

0.16 

0.16 

0.17 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

by (MeV) 

15.8 

15.6 

16.0 

16.0 

17.0 

15.9 

15.9 

15.9 

15.9 

15.9 

15.9 

K (MeV) 

217 

215 

230 

230 

245 

230 

234 

234 

234 

226 

234 

J (MeV) 

30.0 

26.8 

32.0 

32.0 

32.0 

32.0 

28.0 

32.0 

30.0 

30.0 

30.0 

L (MeV) 

47.5 

37.7 

48.3 

45.9 

47.4 

47.2 

7.5 

59.5 

42.0 

35.5 

37.0 

C+ (MeV-fm®) 

57.6 

43.9 

59.3 

60.1 

54.1 

52.7 

49.6 

51.8 

50.9 

51.4 

50.7 

C- (MeV-fm®) 

-4.79 

-0.94 

-22.8 

-24.2 

-15.6 

-13.4 

19.6 

26.0 

36.9 

30.6 

21.9 

Csym 

3.24 

6.07 

1.58 

1.54 

1.77 

1.95 

1.48 

1.40 

1.13 

1.22 

1.46 

P 

-0.64 

-0.54 

-0.58 

-0.52 

-0.62 

-0.65 

-0.48 

-0.47 

-0.51 

-0.48 

-0.48 


TABLE I. Basic parameters of some critical Skyrme forces from Refs. [37, 41], including the L derivatives [20, 41, 44]. In 
addition to these standard quantities, are the isoscalar and isovector constants C± of the energy density gradient terms [Eqs. 
(2) and (8)]; Csym is given by Eq. (12) and the spin-orbit constant j3 is defined below Eq. (9). 



SkM* 

SGII 

SLy5 

SLy5* 

SLy6 

SLy7 

SVsym28 

SVsym32 

SVmas08 

SVK226 

SVkap02 

fcs,o(MeV) 

-2.47 

-0.53 

-12.6 

-13.1 

-9.03 

-7.09 

11.4 

15.6 

37.1 

23.7 

12.7 

^^(MeV) 

-2.48 

-0.46 

-14.6 

-15.0 

-10.1 

-7.61 

13.3 

18.2 

46.7 

29.5 

14.8 

<^0 

163 

21.9 

0.59 

0.92 

1.21 

1.99 

0.90 

0.84 

0.89 

0.79 

0.89 

V 

2.27 

1.89 

0.28 

0.60 

0.62 

0.73 

0.58 

0.61 

0.86 

0.70 

0.59 

Qo(MeV) 

59642 

29908 

73 

72 

137 

287 

-62 

-55 

-62 

-30 

-63 

Q(MeV) 

823 

2570 

42 

41 

63 

98 

-34 

-34 

-34 

-21 

-36 

Toll 

0.006 

0.004 

0.41 

0.43 

0.26 

0.16 

0.43 

0.53 

0.040 

0.89 

0.45 

tII 

0.055 

0.014 

0.59 

0.60 

0.40 

0.28 

0.62 

0.73 

1.68 

1.18 

0.64 

L»o(MeV) i^^Sn 

89 

91 

101 

89 

104 

102 

78 

79 

81 

77 

84 

D(MeV) ®®Ni 

91 

92 

100 

88 

104 

95 

79 

80 

83 

78 

85 

i^^Sn 

89 

91 

100 

89 

103 

95 

77 

78 

81 

76 

83 

208 pb 

90 

91 

109 

88 

102 

93 

77 

78 

81 

76 

82 


TABLE 11. The isovector energy kg and the stiffness Q coefficients are shown for several Skyrme forces ]37, 41, 44]; v is the 
constant of Eq. (27); t/7 is the neutron skin thickness calculated by Eq. (25) with the corresponding L ; the functions D{A) 
for the ELD model in the last three lines are calculated with the relaxation time T having the constant of its frequency 
dependence Tops = 300MeV^- s as explained in the text and in the Figures [51]; the quantities kg ^, Vq, Qo, Tq and Do are 

calculated with L = 0. 
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FIG. 1. (Color online) Isovector W- (13) (with the relevant value of L Refs. [20, 41, 44]) and without (L = 0) derivative L 
constant, and isoscalar w = w+ (see [7]) particle densities are shown vs x = ^ja for the Skyrme force SLy5* (x fa {r — R)/a for 
small nuclear deformations [31, 36, 43]). 
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FIG. 2. (Color online) Isovector density W- (x) (13) (in the logarithmic scale) as function of x within the quadratic approximation 
to e+[e(u))] for several Skyrme forces [20, 37, 41, 43, 44]. 
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FIG. 3. (Color online) IVDR strength functions S(cj) vs the excitation energy hw are shown for vibrations of the nucleus ^®^Sn 
for the Skyrme force SLy5* by dots and dashed lines at L = 50 MeV and solid lines for L = 0; red or magenta (“out-of-phase”), 
and green (“in-phase”) curves show separately the main and satellite excitation modes, respectively (section 4 and 5); the 
collision relaxation time T = 4.3 • 10“^^ s in agreement with the IVGDR widths [51]. 



FIG. 4. (Color online) The same total and different modes (main and satellite) strengths as in Fig. 3 are shown for different 
L = 0 and 60 MeV for the Skyrme force SVsym32. 
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FIG. 5. (Color online) The total IVDR strength functions S{uj) vs the excitation energy huj (in MeV) for different double magic 
nuclei for SLy5* and SVsym32 forces; a slight dependence on the slope parameter L (in MeV) as compared to the L — 0 case 
at the main peaks is shown. 
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(N-Z)/A 


FIG. 6. (Color online) The IVDR splitting vs the asymmetry parameter I = {N — Z)/A. Top: the energies En of the main 
peak (n = 1) and satellite (n = 2); black (n = 1) and red (n = 2) open squares are obtained from the experimental data for 
the integral cross sections of Refs. [54-58] for several Sn isotopes, as explained in the text [30]; the solid and dotted black lines 
with opened circles are the main L=50 MeV and 0 peaks for SLy5* , respectively; the same red curves show the satellites; 
the solid and dashed blue lines with open diamonds denote the main L=60 MeV and 0 peaks for SVsym32* , respectively; 
the same margenta curves show the satellites; the green dots stands for the averaged (IVGDR) experimental data (three last 
from Adrich et al. [32]) and arrows show the particle number of the Sn isotopes. Middle: The ratio of the strengths at the 
satellite to those at the main peak, Sn = S{uJn) (black solids and red dots for SLy5*; blue solids and margenta dashed curves 
for SVsym32); Bottom: The normalized to 100% as explained in the text with the same notations as in the top plot. 
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FIG. 7. (Color online) The IVDR main out-of-phase (dp-, “out”) and in-phase (5p+,“in”) transition densities p-{x) [Eq. (33)] 
multiplied by {r/R)'^, vs x = ^/a « {r — R)/a (spherical nuclei) for the satellite in ^®^Sn with the Skyrme forces SLy5* [43] 
(upper panel) and SVsym32 [41] (lower panel); the two characteristic values L = 0 and L = 50 (or 60) MeV are shown; the 
relaxation time T is the same as in Fig. 3. 










20 



FIG. 8. (Color online) The same but for the IVDR neutron (n) and proton (p) transition densities p^{x) [Eq. (33)] multiplied 
by (r/R)'^, vs x = « (r — i?)/a for the satellite at the energy E 2 in ^®^Sn with the Skyrme forces SLy5* [43, 44] (upper 

panel) and SVsym32 [41] (lower panel); the two characteristic values L = 0 and L = 50 (or 60) MeV are shown too; the 
relaxation time T is the same as in Fig. 3. 
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FIG. 9. (Color online) The IVDR n-p transition densities p“{x) multiplied by {r/R)^ vs the dimensionless distance parameter 
a: = ^/a « (r — R)/a for the same double magic nuclei slightly depending on the slope parameter L (in MeV) for a given 
example SLy5* of the Skyrme forces as compared to the L = Q case near the ES edge as in Fig. 8. 
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FIG. 10. (Color online) The same as in Fig. 9 but for the Skyrme force SVsym32. 
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L (MeV) 


FIG. 11. (Color online) Neutron skin thickness Ar„p = a/s/S (i?„ — Rp) — -^/S/S Tqt ( ro = 1.14 fm) as a function of the 
derivative constant L for the same isotopes as in Table II and Figs. 5, 9 and 10 for the SLy5* and SVsym32 forces; full symbols 
show SLy5* and open ones correspond to SVsym32 calculations; arrows show approximately the values of L for different Skyrme 
forces taken from Refs. [41, 44]. 





